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Method for determining the operating parameters of an 
system comprising a cooling chamber 

The present invention relates to a procedure for 
5 determining operating parameters of an installation for 
thermally cooling articles. 

The invention applies in particular to installations 
for deep- freezing food articles. 

Known deep- freezing installations comprise, for 
example, a deep- freezing chamber or tunnel through 
which a belt conveyor passes, the articles to be frozen 
being deposited on said * conveyor, which runs 
15 continuously or sequentially through the deep- freezing 
tunnel . 

A . cryogenic chamber uses an inert low- temperature fluid 
which exchanges heat directly by contact with the 
20 articles to be deep-frozen. 

Conventionally, a cryogenic chamber uses either dry ice 
(at -80°C) , liquid air or liquid nitrogen (at -196°C) 
as refrigerant. Dry ice allows fresh or deep-frozen 
25 products to be transported without fear of rupture of 
the refrigeration line. Liquid nitrogen and liquid air 
allow either the individual deep- freezing of food 
products, or the hardening of delicate, deformable or 
sticky products (of the ice-cream type, etc.) . 
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The operating parameters of the installations form 
recipes created experimentally. A recipe stores the 
control parameters of an installation for a given 
production run. 

At the present time, there is no sensor capable of 
continuously and contact lessly measuring the internal 
temperature of an article in such a way that the 
recipes can be determined by procedures of determining 
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operating parameters that are then adjusted by means of 
tests performed on site, which tests are often 
destructive . 



5 The procedures of determining recipes generally include 
a step of determining a setpoint for the exit 
temperature of the articles followed by a step of 
determining initial operating parameters and then the 
repetition of a test cycle, comprising a step of 
10 predicting the exit temperature of the articles, a step 
of comparing the predicted temperature with the 
setpoint and, should there be a difference, a step of 
modifying the operating parameters. 

15 The test cycle is repeated until the operating 
parameters result in a predicted temperature 
substantially close to the setpoint temperature. 

It should be noted that all these steps are deduced and 
20 carried out by the operator experimentally, taking into 
account his know-how, his experience and his knowledge. 

If the system consisting of the chamber and the 
mechanism of loading the articles is examined, several 

2 5 parameters may have an influence on the exit 

temperature of the article, namely the production rate, 
which, for a given loading factor, involves a variation 
in the residence time in the chamber, the flow rate of 
the fluid, which acts on the temperature profile, the 

3 0 entry temperature of the articles, the convective 

profile of the chamber, and the loading factor. 

The system is therefore a mult i -variable system and 
existing procedures of determining the parameters, and 
35 especially the step of predicting the exit temperature 
of the articles, do not take these elements into 
account . 

In the prediction steps of the prior art, so as to 
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handle single-variable systems, it has been necessary 
to consider both the convective profile and the loading 
factor as constants and to fix the production rate, the 
temperature of the articles at the inlet of the deep- 
5 freezing chamber, and also other operating parameters 
of the installations. 

Consequently, the recipes determined by the existing 
procedures of determination are relatively imprecise 
10 and require the use of production of articles followed 
by destructive tests thereon. 

The present invention aims to remedy this problem by 
defining a procedure for determining the operating 
15 parameters that is precise and easily implemented. 

For this purpose, the subject of the invention is a 
procedure for determining the operating parameters of 
an installation for the thermal cooling of articles, 
20 comprising a chamber through which said articles pass 
from an inlet to an outlet, using a cooling fluid,, the 
procedure comprising :' 

- a step of determining a temperature setpoint 
for articles at the outlet of said chamber; 

25 - a step of determining initial operating 

parameters for said installation; and 

- a test cycle for testing the operating 
parameters, which comprises: 

- a step of predicting the • temperature of the 
30 articles at the outlet of said chamber, 

- a. step of comparing the setpoint temperature 
with the predicted temperature and 

- if said .comparison step reveals a difference 
greater than a predetermined threshold, a step of 

35 modifying the operating parameters of said 
installation, and a repeat of the test cycle, 
said prediction step being carried out on the basis of 
operating parameters of said chamber, of thermodynamic 
and physical characteristics of said chamber and of 
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thermodynamic and physical characteristics of said 
articles. 



According to other features: 
5 - said prediction step includes a step of 

predicting the behavior of said chamber based on the 
solution of heat balance equations on elementary volume 
slices of said chamber, which is performed at least on 
the basis of thermodynamic characteristics of said 

10 cooling fluid and ' thermodynamic and physical 
characteristics of said chamber; 

said step of predicting the behavior of said 
chamber is furthermore carried out on the basis of 
operating parameters of said installation; 

15 - said operating parameters of said installation 

represent at least one of the elements chosen from the 
group consisting of : 

the speed of a conveyor for transporting said 
article's through said chamber; 

2 0 - the loading factor; and 

the ventilation of the atmosphere in said 

chamber; 

said prediction step includes a step of 
predicting the behavior of said articles based on 
25 solving the discretized heat conservation equation, 
applied to a grid of spatial and temporal points 
constituting a mesh of said articles, which is carried 
out at least on the basis of thermodynamic and physical 
characteristics of said articles; 
30 - said step of predicting the behavior of said 

articles is furthermore carried out on the basis of 
operating parameters of said installation; 

said operating parameters of said installation 
include the temperature of said articles at the outlet 

3 5 of said chamber; 

said step of predicting the behavior of said 
articles is optimized by calculations involving the 
modification of said mesh of said articles using 
mathematical series ; 



said step of predicting the behavior of said 
articles is optimized by omitting the prediction 
calculations for spatial and temporal points of said 
mesh of said articles for which the enthalpy changes 
are below a predetermined threshold value; 

said step of predicting the temperature of said 
articles at the outlet of said chamber is based on said 
step of predicting the behavior of said chamber and on 
said step of predicting the behavior of said articles; 

said step of modifying the operating parameters 
includes a step of manually modifying at least some of 
the operating parameters; 

said step of modifying the operating parameters 
comprises automatically modifying at least some of said 
operating parameters; and 

-* said step of modifying the operating parameters 
comprises modifying at least one of the parameters 
chosen from the group consisting of: 

the flow rate of said cooling fluid; 

the residence time of said articles in said 

chamber; 

the flow rate of gas extracted from ■ said 
chamber; 

the gas speed-up; 

the gas recirculation; and 

the balance between the amounts of incoming air 
and the amounts of outgoing gas. 

- The invention will be more clearly understood on 
reading the description that follows, given solely by 
way of example and with reference to the appended 
drawings in which: 

figure 1 shows a simplified diagram 
illustrating a cooling installation; 

figure 2 is a general flowchart" of the 
procedure for the invention; 

figure 3 illustrates the numerical modeling of 
the articles to be treated; 

figure 4 illustrates the numerical modeling of 



the cooling chamber; and 

figure 5 shows the detailed flowchart of the 
test cycle of the procedure for the invention. 

Figure 1 shows a conventional installation for the 
treatment of food products, the operating parameters of 
which are determined by a procedure according to the 
invention . 

This installation comprises a cryogenic chamber or 
tunnel 2, of conventional type, for freezing food 
articles P by contacting them with a cryogenic fluid 4 
conveyed via a feed line 5 from any source. 

For example, the chamber 2 has a rectangular 
parallelepipedal shape. 

As mentioned above, the cryogenic fluid 4 used may, for 
example, be. dry ice or liquid nitrogen and it is 
injected at one or more points into the chamber 2. 

This chamber 2 is associated with a conveyor 6, of 
conventional type, for introducing the articles P into 
the chamber 2 and for extracting them, said conveyor 
operating either sequentially or continuously. 

The installation has several operating parameters, 
namely the temperature profile in the chamber, the 
residence time of the articles P in the chamber 2 or 
the run speed of the conveyor 6, and the entry 
temperature of the articles P. 

Finally, the installation includes means 12 for 
controlling the quantity of cryogenic fluid 4 injected 
into the chamber 2 . 

These means 12 include means 14 for controlling the 
flow rate of the cryogenic fluid 4. For example, the 
control means 14 consist of systems of solenoid valves 
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or proportional valves of conventional type, placed on 
the feed line 5 for delivering cryogenic fluid 4. 

Advantageously, the installation also includes a gas 
5 ventilation system, controlling, the gas flows, and the 
ventilation of the atmosphere in the chamber 2. 

, For example, this system is composed of specific 

ventilators, for speeding up the gases, of ventilators, 
10 which control the recirculation of the gases, and of a 
combination of ventilators and moving doors, which 
controls the balance between the amount of incoming air 
and the amount of outgoing gas . 

15 The general flowchart of the procedure for determining 
the operating parameters according to the invention 
will now be described with reference to figure 2. 

This procedure starts with a step 1,6 of entering an 
20 exit temperature setpoint for the articles after 
thermal cooling. 

The step 16 is followed by a step 18 of determining the 
initial operating parameters. The parameters determined 

2 5 during this step 18 are known parameters, such as 
mechanical characteristics of the chamber 2 or the 
physical and thermodynamic characteristics of the 
articles P, and variable parameters, such as the 
operating parameters of the installation, which are set 

30 arbitrarily. 

The procedure then includes a step '2 0 of predicting the 
temperature of the articles P at the outlet of the 
chamber 2 . 

35 

This step 20 includes a step 22 of predicting the 
behavior of the chamber 2 and a step 24 of predicting 
the behavior of the articles P. 



Step 22 of predicting the behavior of the chamber 2 is 
used to predict by calculation, as will be described 
later with reference to figure 4, the theoretical 
temperature profile of the cryogenic fluid inside the 
chamber . 2 . 

The results delivered by step 22 depend on the 
thermodynamic characteristics of the cryogenic fluid 4, 
on the convective characteristics of the chamber 2, on 
the characteristics of the means of injecting the 
cryogenic fluid 4, on the characteristics of the 
ventilation system and on the physical characteristics 
of the chamber 2 . 

Step 22 also takes into account the operating 
parameters of the installation, such as the speed of 
the conveyor 6 . 

Step 24 of - predicting the behavior of the articles P is 
used to determine by calculation, as will be described 
later with reference to figure 3, enthalpy changes of 
the articles P according to their external environment 
and to their initial temperature. 

Thus, the results delivered by step 24 of predicting 
the behavior of the articles P depend on their physical 
and thermodynamic characteristics. 

Step 2 2 of predicting the behavior of the chamber and 
step 24 of predicting the behavior of the articles P 
are coupled to each other, as will be described in more 
detail with reference to 5 , so as to deliver a 
theoretical temperature of the articles P at the outlet 
of the chamber 2 . 

Thus, the prediction made in step 20 of predicting the 
temperature of the articles P at the outlet of the 
chamber 2 takes into account the thermodynamic and 
physical characteristics of the chamber 2 and of the 



articles P, and also the operating parameters of the 
installation. 

Consequently, the determination of the temperature of 
the articles P at the outlet of the chamber 2 is 
dynamic, parameterizable and very accurate. 

Step 20 of predicting the temperature of the articles P 
at the outlet of the chamber 2 is followed by a step 26 
of comparing the temperature predicted during step 2 0 
with the temperature setpoint determined during 
step 16. 

For example, this comparison step 2 6 takes into account 
a tolerance interval of the order of a few degrees 
around the exit temperature setpoint determined during 
step 16. 

When the difference between the predicted temperature 
and the setpoint is above a predetermined threshold 
value, comparison step 26 is followed by a step 28 of 
modifying the operating parameters of the installation. 

For example, the operating parameters modified during 
this step 2 8 include the speed of the conveyor 6 and 
the parameters involved in determining the theoretical 
temperature, profile of the fluid 4 in the chamber 2, 
i.e. for example the flow rate of the cryogenic fluid 
4, the control of the ventilation inside the chamber 2 
and the loading factor of the conveyor 6. 

The modifications of the operating parameters during 
step 28 may be carried out directly by an operator or 
automatically by a computer, taking the maximum and 
minimum limits for each of the parameters into account. 

In addition, the modifications may each time affect one 
or more parameters . 
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It is also possible to define an order of modification 
of the operating parameters in order to try to achieve 
the setpoint temperature by modifying only one 
parameter at a time. If the setpoint cannot be achieved 
5 by modifying a first parameter between limit values, 
this parameter is set to a limit value or a mean value, 
and in the following iterations the next parameter in 
the list is modified. 

10 After step 28. of modifying the operating parameters, 
the procedure returns to step 20 of predicting the 
temperature of the articles P at the outlet of the 
chamber 2, the procedure thus forming a test cycle for 
testing the operating parameters of the installation, 

15 comprising step 20 of predicting the exit temperature 
of the articles, step 26 of comparing the predicted 
temperature with the setpoint temperature and step 2 8 
of modifying the operating parameters. 

2 0 This test cycle, denoted by the general numerical 

reference 29, is repeated until the comparison made 
during step 2 6 . between the predicted temperature and 
the setpoint temperature reveals a difference below a 
predetermined threshold value. 

25 

The cycle 2 9 is then interrupted and step 2 6 is 
followed by a step 3 0 of recording the last operating 
parameters tested, which thus form a recipe. 

3 0 Since the prediction made during step 2 0 is very 

accurate and takes into account the operating 
parameters of the chamber 2, the thermodynamic and 
physical characteristics of the chamber 2 and the 
thermodynamic and physical characteristics of the 
35 article P, the recipe determined by the procedure for 
the invention is accurate and close to the actual 
operation . 

In addition, such a recipe is easily adapted to 



modifications in the operating conditions. 

For example, if an operating parameter such as the 
temperature of the articles on entering the chamber 2 
is modified, the recipe may be corrected by executing 
the procedure for the invention using the actual recipe 
parameters as initial operating parameters during step 
18, the execution of the procedure for the invention 
making it possible for the corrections to be made to 
the operating parameters, in order to obtain the proper 
exit temperature of the articles P, to be rapidly and 
simply determined. 

Step 20 of predicting the temperature of the articles P 
at the outlet of the chamber 2 will now be described in 
greater detail with reference to figures 3, 4 and 5. 

Figure 3 shows an example of the mesh for a food 
article P as employed during step . 24 of predicting the 
behavior of the articles P. 

In the procedure for determining the operating 
parameters, the thermodynamic and physical 

characteristics of the articles P are taken into 
account in step 24 of predicting their behavior on the 
basis of a modeling of the articles P to which the 
discretized heat conservation equation is applied. 

This is because the heat conservation equation cannot 
be solved at any point in space and at any instant by a 
simple integral function. 

The procedure employed consists in discretizing this 
equation in such a way that it is solved only at the 
spatial and temporal points called nodes, these being 
denoted by the general numerical reference 32. 



After defining a mesh for the article P, the heat 
conservation equation is applied to each of the 
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nodes 3 2 . 



What is thus obtained is a system of equations that 
have to be solved in order to determine the thermal 
state of the article P both in time and in space: 



dxy dx J ay\ oy j 



dz p dt 



in which X, Y and Z are axes that define an orthonormal 
10 spatial reference frame around the article P, T is the 
temperature of the article P expressed in kelvin (K) , 
and C is its specific heat expressed in watts per 
kilogram per kelvin (W/(kg.K)) . 

15 The food articles P to be deep-frozen generally consist 
of different bodies. 

This means that the phase change is accompanied by a 
change in temperature, and that the .heat conservation 
2 0 equation may still apply. 

However, if it is necessary to treat a pure substance, 
the equation is no longer continuous. In this case, the 
problem is simplified by modifying the enthalpy table 

2 5 of the pure substance so that the phase change 

generates a small change in temperature. 

The discretization is carried out using the 
mathematical method of finite differences in variable 

3 0 mode . 



As is known, this may be carried out in two ways. 

The first way - implicit discretization - has the 
35 advantage of being stable whatever the spatial and 
temporal configuration. At a given instant, it allows 
the temperature of a node 32 to be determined as a 
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function of the temperature of the neighboring nodes at 
the same instant. However, it does involve constant 
boundary conditions and a matrix solution of the system 
of equations formed by each of the nodes 32. 

5 

The second way - explicit discretization - allows 
direct determination of the temperature of a node 3 2 at 
an instant T+AT on the conditions at the instant T. The 
result is immediate, however it is necessary to choose 
10 a suitable step time so as to avoid any instability of 
the model . 

The first method is recommended if it is desired to 
obtain mainly the surface temperature of an article, 
15 which corresponds to the operation commonly referred to 
as "crust freezing" . The second method is recommended 
when it is desired to deep- freeze an article and to 
know the temperature in the core thereof. 

20 The mesh for the article P is a crucial problem. On it 
will depend directly the simplicity of the subsequent 
processing and the precision of the results. 

A large number of nodes that provide great precision in 
25 the result, but requires a long computation time. A 
compromise has to be found between precision and 
computing time. 

For example in the case of a food product having 
30 external dimensions of 100 x 60 x 10 mm with a regular 
mesh of cells every millimeter, more than 17000 nodes 
are required and the same number of equations in order 
to define the behavior of the article P. 

35 In the described version of the invention, step 24 
includes calculations that make it possible to optimize 
the mesh used for predicting the behavior of the 
products P. 
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For example, in the case of crust freezing it is more 
particularly solidification of thin skin of the product 
caused by a phase change that is monitored. It is 
therefore necessary to have a dense mesh around the 
periphery, but a coarser one in the core . 

In order not to manually input the coordinates of each 
of the nodes and to maintain simple relationships 
between the nodes and to facilitate the processing, one 
solution consists in distributing the nodes in each 
direction in space using, for example, a geometric 
progression, as shown in figure 3 . 

For example, along the X axis, the nodes are 
distributed as follows: taking Ax as the value of the 
first term, which corresponds to the abscissa of the 
first node, and r as the common ratio, which differs 
from 1, of the geometric series used, then the value of 
the nth term is Axr 11 " 1 , which corresponds to the 
position on the X axis of the nth node. The sum of the 
n first terms is: 

i- r " 

S = Ax + — 

i-r " 

Figure 3 shows the positions of the nodes according to 
this mesh on a parallelepipedal article P on which a 
parity condition on the number of nodes has been 
imposed so as to simplify the solution. 



The value corresponding to the length of the product P 
along the X axis is thus obtained: 



- 15 - 



where R = 



1 



1 + ; 



with 1, which corresponds to the abscissa of the 
central node on this length, given by: 



f f 
1- 
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u 



Ax 



l4t ) 



The imprecision along^ the X axis is then expressed as 
follows : 

10 



-1/2 < /- 




To result in simple calculations, the imprecisions on 
the three axes are set to the same value. This 
15 introduces an error in the dimensions of the product P 
that is acceptable if one is only interested in the 
temperatures over a thin skin and if the core 
temperatures vary little, as will be the case in crust 
freezing operations . 

20 

In the case of deep- freezing operations, in which the 
aim is to determine the core temperature of the 
product; a corrective term may be inserted into the 
formulae. In the case of the X axis, the following 
25 corrective term is inserted: 
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Another possible optimization method consists in 
reducing the processing time by omitting certain 
5 calculations. 

Specifically, on each node, the heat flux is added to 
the six faces of its elementary volume. However, there 
are zones in which the thermal effects can be likened 
10 to one-dimensional problems. 

.To exploit this particular feature, the processing is 
broken down by summing no longer the overall thermal 
fluxes on each face, but in each direction. A time 
15 interval AT, the equations at the nodes are solved in 
each direction by going from the boundary towards the 
core, until the enthalpy change is considered as being 
negligible, since it is below a predetermined threshold 
value . 

20 

By carrying out this operation in each direction, a 
volume of the product P encompassing all the nodes for 
which the enthalpy changes will be negligible, and 
therefore for which no calculation will be carried out, 
25 is defined. 

It is thus possible to save computing time, in 
particular over the first few moments of the exchange. 

30 If the article P is of complex shape, It may be 
decomposed into a set of elementary shapes to which the 
mesh defined above, or any other mesh suitable for the 
shape of the article P, is applied. 

35 Figure 4 shows schematically the chamber 2 for treating 
the food articles P, as modeled for implementing step 
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22 of predicting its behavior. 

In the procedure for determining the operating 
conditions, the thermodynamic and physical 

5 characteristics of the. chamber 2 are taken into account 
during step 22 of predicting the behavior of the 
chamber 2 on the basis of its model in the .form of 
elementary slices . 

10 As described previously with reference to figure 1, the 
cooling chamber 2 is associated with a conveyor 6. The 
chamber is supplied with cryogenic fluid 4 via a feed 
line 5. The chamber 2 can be likened to a rectangular 
parallelepiped . 

15 

To determine the theoretical temperature profile of the 
fluid 4, the method used to predict the behavior of the 
chamber 2 consists in carrying out a succession of 
local heat balances. 

20 

For this purpose, the thermodynamic system of the 
chamber' 2, in the steady state, is modeled in the form 
of elementary slices 34 x to 34 n perpendicular to its 
length. The sum of these elementary slices 34 x to 34 n 
25 represents the internal volume of the chamber 2. 

For each elementary slice 34 x to 34 n/ the heat transfer 
balance is made, so as to determine the enthalpy of the 
fluid 4 and therefore its temperature. 

30 

This balance must take into account the following: 

- the thermal losses with the outside of the chamber 2; 

- the cryogenic liquid 4 injected into the spray zones; 
and 

3 5 - the heat exchange between the products P and the 
fluid 4 . 



In the case of the slice 34i of the chamber 2, having 
dimensions Lxlxh, the thermal balance is given by the 
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following equation : 




H jiO) H MO 



In this equation: 

H fS (i) corresponds to the enthalpy of the cryogenic 
fluid 4 on leaving the elementary slice 34 i# expressed 
in joules per kilogram (J/kg) ; 

H fed) image of 14/1] corresponds to the enthalpy of 
the cryogenic fluid 4 on entering the elementary slice 
34 i7 expressed in joules per kilogram (J/kg) ; 
H fLiq corresponds to the enthalpy liquid of the 
injected cryogenic fluid 4, expressed in joules per 
kilogram (J/kg) ; 

H P e(ij corresponds to the enthalpy of the article P on 
entering the slice 34 i# expressed in joules per 
kilogram (J/kg) ; 

H 'p S (i) corresponds to the enthalpy of the article P on 
leaving the slice 34±, expressed in joules per kilogram 



K T corresponds to the heat exchange coefficient for 
exchange between the chamber 2 and the outside, 
expressed in watts per square meter per kelvin 
(W/ (m 2 K) ) ; 

TWys^qKi) corresponds to the mass flow rate of the 

cryogenic fluid 4 vaporized in the slice 34 ±/ expressed 
in kilograms per second (kg/s); 

TJ'lfain corres P onds to the mass f l° w rate of cryogenic 
fluid 4 entering the slice 34 i7 expressed in kilograms 
per second (kg/s); 

772 p corresponds to the mass flow rate of products to be 

treated, expressed in kilograms per second (kg/s); 

T Amb corresponds to the ambient temperature, expressed 

in kelvin; and 

T f e<i) corresponds to the temperature of the cryogenic 



(J/kg) ; 
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fluid 4 entering the slice 341, expressed in kelvin. 

Figure 4 also shows the following thermal fluxes: 
mjSproytfiH jz iq is represented by the letter A; 

5 thfe® H Mn is represented by the letter B; 
fhjsftll mo is represented by the letter C; 
lh P Hpe(i) is represented by the letter D; and 
lllrpH. ps(J) is represented by the letter E, 
where 171 j&(o"~/W Mii^^W J5p*-ay{i) . 

10 

From experience, it is known that, under certain 
operating conditions (too low a production rate or too 
low a temperature of the cryogenic fluid 4) , the 
injected cryogenic fluid 4 is only partly vaporized and 
15 a fraction of the liquid flows away toward the inlet of 
the chamber 2 . 

If it is desired to take this phenomenon into account, 
it is preferable to solve the local balance equations 
20 by starting with the elementary slice located at the 
outlet of the chamber 2. The calculations are therefore 
performed in the opposite sense to the path of the 
articles P, i.e. along the X axis as shown in figure 4. 

25 This is because, in this direction, the fraction of 
unvaporized liquid can be transferred into the next 
slice and so on, until the ventilation zones are 
reached, where the injected flow rates are zero and 
where the surplus liquids are vaporized. 

30 

To determine the fraction of unvaporized cryogenic 
liquid 4 in an elementary slice, a limiting fluid 
enthalpy is denoted, below which it will appear as 
liquid. 
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This amounts to setting a minimum gaseous fluid 
temperature in the chamber 2 . 

The unvaporized liquid content on leaving the 
elementary slice 34i corresponds to x Mi) an< 3 ^ s 
expressed in the following form: 



If the calculations are simplified by considering that 
the enthalpy of this liquid fraction is approximately 
equal to the enthalpy of the injected cryogenic fluid 
4, the following expression for the amount of liquid is 
obtained: 



In this equation, H fL i m corresponds to the limiting 
enthalpy of formation of an amount of liquid in an 
elementary slice of the chamber 2. 

Figure 5 shows the details of the test cycle 29 and in 
particular step 20 of predicting the temperature of the 
articles P at the outlet of the chamber 2. 

To be able to predict the temperature of the articles P 
of the outlet .of the chamber. 2, the cooling procedure 
involves step 22 of predicting the behavior of the 
chamber 2 and step 24 of predicting the behavior of the 
articles P. 
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The test starts by carrying out step 22 of predicting 
the behavior of the chamber 2, during a step 40. 
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This step 40 delivers the heat losses 42 per elementary 
slice, these being reintroduced into step 22. 



After this operation has been reiterated a certain 
5 number of times, the total heat losses 44 and the 
temperature profile 46 of the fluid 4 in the chamber 2 
are obtained during step 40. 

To calculate the heat balance for each slice, step 22 
, 10 requires the enthalpy changes of the articles P. 

Because, during the first iteration, the temperature 
profile of the fluid 4 in the chamber 2 cannot be 
calculated, it is set arbitrarily. 

15 Next, step 24 of predicting the behavior of the 
articles P is carried out during a step 50. This step 
50 delivers the enthalpy 52 of the article P at the 
outlet of the chamber 2, i.e. its temperature. 

20 Optionally, step 24 of predicting the behavior of the 
articles P also delivers the enthalpy changes 54 of an 
article P for each elementary slice of the chamber 2. 
In this case, this datum is sent back to step 22 of 
predicting the behavior of the chamber 2, which inserts 

2 5 it into the thermal balance equation for each 

elementary slice. 

The enthalpy 52 of the article P at the outlet of the 
chamber 2, the temperature profile 46 of the fluid 4 in 

3 0 the chamber 2 and the total heat losses 44 are put into 

the relationship, at step 60, so as to determine the 
total flow rate of the fluid. 

Optionally, the flow rate 62 of fluid injected into 
35 each elementary slice is also obtained. In this case, 
this datum is sent back to step 22 of predicting the 
behavior of the chamber 2, which inserts it into the 
heat balance equation for each elementary slice. 
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Next, a check is made at step 70 to determine whether 
the temperature profile of the fluid 4 in the chamber 2 
is stable. 

5 • For example, the temperature profile of the fluid is 
considered as stable if it meets, twice in succession, 
the following criterion: 



10 



H*o* & P* l > k "» <dif ^profile 



In this equation, dif_jprofile is a constant set by the 
operator . 

The profile is considered as unstable for the first 
15 pass. 

As long as the profile is considered as unstable, the 
procedure returns to step 4 0 and the succession of 
operations for defining a profile is repeated. 

20 

Once a stable profile has been obtained, a check is 
made during step 26 to determine whether the setpoint 
determined during, step 26, relating to the temperature 
of the articles P at the outlet of the chamber 2, has 

2 5 been reached. 

If the setpoint has been reached, the operating 
parameters resulting in the final temperature profile 
of the fluid 4 inside the chamber 2 are recorded during 

3 0 step 3 0 and these form a recipe. 

If the setpoint has not' been reached, the operating 
parameters of the installation are modified at step 28. 
For example, this modification includes a correction 
35 102 to the flow rate of the fluid 4 before reiteration 
of the algorithm. Optionally, the modification includes 



a correction 104 made directly to the operating 
parameters that determine the temperature profile of 
the fluid 4, which is used in step 22 to predict 
behavior of the chamber 2 . 

In this example, the predicted temperature of the 
articles at the outlet of the chamber is used to 
determine the flow rate of the fluid 4 injected into a 
cryogenic chamber 2 . 

In the same way, it is possible to vary the residence 
time of the articles P in the chamber 2 by modifying 
the speed of the conveyor 6 or the stop times in the 
case of a sequential conveyor. It is also possible to 
vary the gas extraction rate or the loading factor, or 
a combination of these parameters. 

Likewise, it is possible to vary the balance between 
the amount of incoming air^ and the amount of outgoing 
gas, the gas extraction rate, the gas speed-up," or the 
gas recirculation by acting on the elements for 
controlling these parameters. 

Moreover, if the procedure of the invention is used for 
an installation employing contactless sensors for 
measuring the exit temperature of the articles, for 
example sensors based on thermal radiation or infrared 
imaging, or else based on an MWT (microwave 
thermometry) measurement, such as the sensor described 
in patent FR-A-2 771 552, the results delivered during 
the step of predicting the temperature of the articles 
at the outlet of the chamber can then be crosschecked 
with the measurements delivered by these sensors. 

In this case, one or other of the information items is 
used to verify the other. 

In one embodiment, the information delivered by the 
sensor is used to correct the prediction. 
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Although one particular embodiment has been described, 
this should not be considered as limiting the scope of 
the present invention. 

5 

As a variant, the cooling procedure of the invention 
may also be applied in a mechanical refrigeration 
installation having an indirect heat exchange device. 

10 The invention has been described in the case of cooling 
food articles, however it may also be applied to other 
types of articles, especially metal articles. 

In addition, the term "cooling" also covers systems 
15 aimed at maintaining and controlling a temperature 
below the initial temperature of an article. 

The procedure of the invention may . be implemented 
using, for example, a program executed on a computer or 
2 0 any other suitable software and/or hardware approach. 



